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(54) Adaptive tonal control system witfi constrained output and adaptation 

(57) An adaptive control system and method for 
actively canceling tones in an active acoustic attenuation 
system has an adaptive parameter bank. Adaptation of 
the adaptive parameter bank can be constrained with 
respect to the null space of a C model of an auxiliary path 
(e.g. speaker-error). Alternatively, output from the adap- 
tive parameter bank can be constrained with respect to 
the effective null space of the C model. The preferred 
system uses singular value decomposition, normaliza- 
tion, and demodulation to implement the methods of con- 
straining output and adaptation. The invention can 
eliminate stability and tracking problems associated with 
over-parameterization. 
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Description 

Field of the Invention 

5 The invention relates generally to adaptive control systems and methods. The invention is particularly useful for 
actively canceling tones in an active acoustic attenuation system. 

Background of the Invention 

10 In many active control applications, cancellation is required only at discrete frequencies where tonal disturt)ances exist. 
While adaptive filtert)ased systems using the filtered-X and fittered-U LMS update methods are effective in canceling 
broadband disturt>ances, adaptive filters can be computationally burdensome when canceling tones, especially if the 
number of cancellation actuators and error sensors is large. 

Moreover, in tonal applications, broadband adaptive filters can sometimes be over-parameterized. Over-parameteriza- 
15 Won can lead to lack of persistent excitation withi n the f liter. Another problem with overparameterizati on is that the potential 
exists for output from separate cancellation actuators to increase yet cancel one another, thus avoiding detection by 
error sensors. This can lead to unnecessarily high power consumption, and instability. These problems can be remedied 
via leakage methods, but most leakage methods can compromise performance. 

Another characteristic of adaptive filters is tfiat, with respect to adaptation, there exists an interdependence among 
20 all tap weight values In the filter. This interdependence can reduce convergence rate. 

As an alternative to the filtered-X and fiitered-U LMS updates. Fast Fourier Transforms can be used to transform 
signals from error sensors into tiie frequency domain as a set of complex numbers. The real and imaginary part can 
then be separately filtered through the complex transpose (i.e. Hermitian transpose) of a transfer function representing 
the speaker-error path. This procedure accounts for phase shifts and delays through the speaker-error patii and can be 
25 used to improve stability But, this procedure can be quite burdensome computationally, and also increases time to track 
changes in tiie system. 

It is therefore desirable to provide a method of attenuating selected tones in a system input that avoids over-param- 
eterization, accounts for phase shifts or delays through the speaker-en-or path, reduces the required amount of process- 
ing and improves the convergence rate of the system, while at tiie same time maintains long term stability. 

30 

Summary of the Invention 

The invention provides an adaptive control system and method that is capable of effectively attenuating selected 
tones in a system input without the problems described above. 

35 In one aspect, the invention is a tonal control system and method in which the output from an adaptive controller is 
constrained with respect to the null space (or an effective null space) of an auxiliary path (e.g. a speaker-error path in 
a sound cancellation system). The system has a plurality of actuators each receiving a correction signal yn and outputting 
a secondary input which combines with a system input to yield a system output. A plurality of error sensors sense tiie 
system output, and each error sensor generates an error signal. An adaptive controller outputs the correction signals 

40 y^ to the plurality of actuators. 

The adaptive controller has an adaptive parameter bank that outputs a plurality of output signals in accordance with 
adaptive parameters. The adaptive parameters in tiie adaptive parameter bank are normally updated using the error 
signals or some signal derived therefrom. The adaptive controller also has a model of the auxiliary path between the 
output of tiie adaptive controller and the error sensors which is referred to as a C model. The C model has an effective 

45 range and an effective null space, the effective null space being defined as tiie subspace spanned by yp, such that Cyp 
= 0 or close to zero where C is a matrix representing the C model, and y^ is a vector representing tiie correction signals. 
The adaptive confroller also has an output weighting element that inputs the output signals from tiie adaptive parameter 
bank and weights the output signals to generate correction signals yp which are constrained to be within the effective 
range of the C model. 

50 A preferred method to weight the output signals so tiiat the correction signals are within tiie effective range of tiie 
C model is to use a singular value deconposition of a C model matrix such that: 

c = usv" 

55 where S is a real-valued matrix in which the off diagonal elements are zero, U is a unitary matrix, and is the Hermitian 
transpose of a unitary matrix V. The output weighting element can then include a matrix representing -VS^U", and it 
may be prefen-ed in many applications to selectively normalize tiie diagonal elements in 8^. In another preferred embod- 
iment, the output weighting element can include a mati-ix representing V. In tiiis embodiment, it is preferred that the 
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adaptive controller also have an error weighting element that includes a matrix representing -S^U^, where again it may 
be preferred to selectively normalize the diagonal elements of S^. 

In another aspect of the invention, the adaptation of the parameter bank is constrained with respect to the effective 
null space of the C model, in such a system, the adaptive controller has an adaptive parameter bank that outputs a 

5 plurality of output signals in accordance with adaptive parameters. The output signals are used to generate the correction 
signals yn, and in embodiments without an output weighting element, the output signals can be used directly as the 
correction signals yn- The adaptive controller has an error weighting element that inputs the error signals from the error 
sensors and weights the en'or signals to generate error input signals used to update the adaptive parameters in the 
adaptive parameter bank. The error weighting element generator en-or input signals such that adaptation Is constrained 

10 to be within the effective range of the C nxxdel. In one preferred system, the error weighting element can be represented 
by the matrix -VS^U^. In another embodiment described above, the output weighting element can be set to V, and the 
error weighting element set to -S^U" where it may be preferred to selectively normalize S^. 

In another aspect, the Invention provides a system and method of demodulation for determining the adaptive param- 
eters in the adaptive parameter bank so that the system operates to minimize a set of error signals. In this aspect, the 

15 adaptive controller in the preferred system has an adaptive parameter bank tiiat inputs a reference signal during sequen- 
tial sample periods and outputs a plurality of modulated output signals in accordance with adaptive parameters for each 
of the sequential sample periods. An error weighting element receives enror signals from en'or sensors and outputs a 
plurality of error input signals. The adaptive controller also has a parameter update generator that applies both an in- 
phase demodulation signal and an quadrature demodulation signal to the error input signals for each sample period 

20 when generating update signals, and outputs an in-phase update signal and a quadrature update signal that are used 
to update the adaptive parameters In the adaptive parameter bank. The in-phase demodulation signal preferably has 
the same frequency as the tone being controlled, and the quadrature demodulation signal is preferably shifted 90^. In 
such a system, the error or the output weighting elements can be -C^ or -C"\ but it Is preferred to constrain output and/or 
adaptation as discussed above. 

25 An object of the invention is to improve system stability and other problems associated witii over parameterization. 
The Invention can accomplish these objects by constraining adaptation of the adaptive parameter bank with respect to 
the effective null space of the C model of the speaker-error path. Alternatively, the Invention can accomplish these 
objectives by constraining the output from the adaptive parameter bank with respect to tiie effective null space of tiie C 
model. 

30 Anotiier object of the invention is to account for the effects of propagation delay through tiie speaker error path. The 
invention accomplishes this object through demodulation which accounts for the effect of any change in actuator output 
to the received error signal. 

Yet another object of the invention is to improve frequency tracking. The invention can accomplish this object by 
matching output from the adaptive parameter bank in consecutive sample periods, even when the output weighting 
35 matrix is replaced, to promote fast tracking. 

Brief Description of the Drawings 

FIG. 1 is a schematic illustration of an active acoustic attenuation system that attenuates a tone at a discrete fre- 
40 quency in accordance witii tiie invention. 

FIG. 2 is a schematic illustration of an active tonal attenuation system that has an error weighting element in accord- 
ance with the invention. 

FIG. 3 Is a schematic illustration similar to FiG. 2 showing another embodiment of an active tonal attenuation system 
having an error weighting element, in accordance with tfie Invention. 
45 FIG. 4 is a schematic illustration of an active tonal attenuation system tiiat has both an error weighting element and 
an output weighting element in accordance with the Invention. 

FIG. 5 is a schematic illustration of an active tonal attenuation system that has an output weighting element in 
accordance with the Invention. 

50 Description of the Preferred Embodiments 

FIG. 1 1llustrates an active acoustic attenuation system designated generally as 1 0. The system 1 0 uses an adaptive 
controller 12 to attenuate a tone at a particular frequency in a disturt^ance 18. The adaptive controller 12 is preferably 
embodied within a programmable digital signal processor. The adaptive controller 12 has an adaptive parameter bank 
55 13, a parameter update generator 28; and erther an output weighting element 14. an error weighting element 26, or both 
14 and 26. To attenuate several tones at distinct frequencies, several attenuation systems 10 such as shown in FIGS. 
1-5 can be Implemented separately and contemporaneously on the same digital signal processor. Separate tones are 
substantially orthogonal so an adaptive confroller 12 implementing separate and contemporaneous tonal attenuation 
systems 10 can effectively attenuate several tones in a disturbance 18. 
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In the adaptive controller 12. the adaptive parameter bank 13 generates a plurality of m output signals y. The m 
output signals y can be applied to an output weighting element 1 4 to generate n correction signals yn. It is preferred that 
the m output signals y be a vector of digital signals, and that the output weighting element 1 4 be an m x n output matrix. 

Each of the n correction signals yn drives an actuator 16 that provides a secondary input or cancellation signal 1 7 
5 that combines with a system input to yield a system output 21 . That is, the secondary inputs 1 7 from the actuators 1 6 
propagate into the system and attenuate the disturbance 1 8 to yield the system output 21 as represented schematically 
by summing junction 20. A plurality of p en-or sensors 22 senses the system output 21 , and generates p error signals 
Cp. In FIG. 1, the path of the n correction signals yp through the n actuators 16, the path of the secondary inputs or 
cancellation signals between the actuators 16 and the error sensors 22. and the path through the p error sensors 22 Is 
10 defined as a p x n auxiliary path (e.g. a p x n speaker-error path), and is illustrated by block 24. 

The adaptive controller 12 receives an error signal Cp from each of the p error sensors 22. The controller 12 can 
have an error weighting element 26 that processes the p error signals e^ to yield m error input signals e. The error 
weighting element 26 is preferably an m x p matrix. 

The parameter update generator 28 in the controller 12 receives the m error input signals e. and generates a set of 
15 parameter updates u. The parameter updates u are used to adapt one or more scaling vectors in the adaptive parameter 
bank 13. A scaling vector can be adapted by accumulating the updates u with the existing scaling vector. The scaling 
vector is then typically applied to a tonal reference signal to generate the m adaptive output signals y. 

The output weighting element 14 and the error weighting element 26 can be chosen to constrain the output from 
the controller 12 (i.e.. constrain the correction signals transmitted to the actuators 16), and/or to improve the conver- 
20 gence of the adaptation process. There are several methods for generating the output 14 and error weighting elements 
26, but it is preferred that a model of the speaker-en'or path 24, referred to herein as a C model, be used to generate 
the output 14 and error weighting elements 26. The C model can be generated off-line, but it is preferred that the C 
model be adaptively generated on-line as described in U.S. Patent No. 4,677.676 which Is incorporated herein by ref- 
erence for the purposes of adaptive on-line C modeling. In the prefen-ed system 10, the C model is represented by a p 
25 x n matrix C. 

The error sensors 22 preferably generate en-or signals ep every sampi e period k. It is desirable to adapt the controller 
12 rapidly In real time with respect to sample period k. This can be approximated over time by den^ulating the error 
input signals e by the in-phase and quadrature components of the particular frequency being attenuated. The demod- 
ulation is accomplished using in-phase and quadrature demodulation signals in the parameter update generator 28. The 

30 in-phase and quadrature components are formed for the particular frequency being attenuated. Therefore, the need for 
performing Fast Fourier Transforms on the error signals has been alleviated, as well as other off-line analysis such as 
averaging or integrating over one or more periods to account br propagation delay in the auxiliary C path. 

FIG. 2 illustrates a system 10 implementing the demodulation method described generally above. In FIG. 2, the 
system 10 implements a weighted error method of adaptation. The system 10 in FIG. 2 is similar to tiie general system 

35 10 shown in FIG. 1 in many respects, expect the output weighting element 14 Is omitted in FIG. 2 (I.e. set to Identity). 
An advantage of omitting the output weighting element 14 is that there is a savings in processing requirements. Also, 
when tracking over large ranges of frequencies it may be useful to have multiple versions of the output 14 and the error 
weighting element 26 available on the digital signal processor, and omitting the output weighting element 14 reduces 
the burden of switching between multiple versions. 

40 In FIG. 2, the adaptive parameter bank 13 In the controller 12 receives an input signal x(k) from an input sensor 30. 
The input signal Is fransmitted to a phase locked loop circuit 32 in the controller 12. The phase locked loop circuit 32 
outputs a reference signal at a particular frequency which is the frequency of tiie tone being attenuated. In particular, 
the reference signal is preferably a discrete time sequence in the form of a cosine wave at a particular frequency. It is 
preferred that tiie reference signal have a normalized magnitude. Other metiiods of obtaining a reference signal can be 

45 used within the spirit of the invention, however, the phase locked loop circuit 32 is preferred because it allows frequency 
tracking. 

The reference signal is separated into two signals at junction 34: An in-phase reference signal is transmitted through 
line 36, and a quadrature reference signal is transmitted through line 38. The in-phase reference signal is ti-ansmitted 
through line 36 to an in-phase scaling element 40. The in-phase scaling element 40 multiples the In-phase reference 

50 signal by an in-phase scaling vector Yr to generate m in-phase components y^ of tiie adaptive output signals yp. The 
in-phase scaling element 40 stores tiie values of tiie in-phase scaling vector Yr . and updates tiie values. The values 
of Yr are updated by summing the product of an in-phase update signal Ur multiplied by a step size ^. 

Contemporaneously, quadrature components yj of the output signals y„ are generated. The quadrature reference 
signal is transmitted through line 38 to a phase shifter 42 that shifts the quadrature reference signal 90** to In effect 

55 generate a sine wave coaesponding to the cosine wave. Thus, in this context, the term quadrature reference signal 
corresponds to a reference signal that has been phase shifted 90" from the in-phase reference signal. The quadrature 
scaling element 44 multiplies the quadrature reference signal by a quadrature scaling vector Y| to generate m quadrature 
components yj of the adaptive output signals yp. The scaling element 44 stores the values of the quadrature scaling 
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vector Y|. and updates the values by summing the values by the product of a quadrature update signal Uj multiplied by 
the step size 

The m in-phase output signals yr and the m quadrature yj output signals are summed at summer 46 to generate m 
con-ection signals yn. where in the case of FIG. 2 m » n . The n con-ection signals yn are transmitted to n actuators 16. 

The an-ay of error sensors 22 generate p error signals e^. The p en'or signals Cp are transmitted to error weighting 
element 26. The en'or weighting element 26 can be determined using the p x n C matrix to eliminate problems associated 
with over-parameterization and to also account for phase shifts and delay in the auxiliary C path. The matrix C can be 
decomposed at the frequencies of interest using singular value decomposition as represented t^elow: 



where U is an p x p matrix, S is a p x n matrix, and is an n x n Hermitian transpose of an n x n matrix V. The matrices 
U and V are unitary matrices, and the off diagonal elements of S are zero while the diagonal elements are in general 
real and positive. In one embodiment, the error weighting element 26 preferably applies the n x p matrix H g = -VN ^ U ^ , 
where V is an n x n matrix; is an n x p matrix that is the Hermitian transpose of normalizing matrix N which is formed 
by inverting some of the values on the diagonal of S (e.g.. the values that are not zero or dose to zero): and is the 
Hermitian transpose of an p x p matrix U. 

The C matrix has an effective range and an effective null space. The effective range is the subspace spanned by 
the columns of V corresponding to relatively large singular values (e.g. all singular values larger than .02 times the largest 
singular value). The effective null space con-esponds to the subspace spanned by the remaining columns in V (e.g. the 
smallest singular values). In other words, the effectve null space is defined by: 



where C Is a matrix representing the C path model and Yp is a vector representing non-trivial correction signals. 

One purpose of applying -VN^U^ is that adaptation will not occur in the effective null space of C. Since the columns 
of V corresponding to the effective null values of N are not included in the adaptation process, any components that 
previously existed or that accumulate as a result of noise in the system are likely to be left unchecked. Therefore, it is 
prefen*ed in this emtxxJiment to subtract or leak from the output signals yr and yj as needed (i.e. conrtponents of y In the 
effective null space of C). 

The use of normalizing matrix N can improve the rate of convergence. In the case of n = p [ and N = S"\ the n x 
pmatrix = -VN^U^ is equal to -C""'. Alternatively, the errorweighting element 26 can apply an nxp matrix H 2 = -C^ 
where is the Hermitian transpose of the p x n C matrix. 

Error weighting element 26 preferably has a junction 48, an in-phase weighting element 50 and a quadrature weight- 
ing element 52. Each of the p error signals ep is transmitted to the junction 48, and the p error signals ep are then 
contemporaneously transmitted to the in-phase weighting element 50 and to the quadrature weighting element 52. The 
in-phase element 50 of the error weighting element 26 contains the real parts of the complex elements of the error 
weighting matrix H2. The quadrature element 50 of the error weighting element 26 contains the coefficients of the imag- 
inary parts of the complex elements of the error weighting matrix H2. Both the in-phase 50 and the quadrature 52 elements 
of the error weighting element 26 contain real values. When referring herein to in-phase and quadrature weighting ele- 
ment, the term in-phase weighting element refers to the real parts of tiie complex elements in a weighting matrix, and 
the term quadrature weighting element refers to the imaginary parts of the complex elements in a weighting mat-ix. The 
p en-or signals Cp are processed contemporaneously through the in-phase element 50 and the quadrature element 52 
to each provide m error input signals e (i.e. n = m in the case of FIGS. 2 and 3). Botii sets of m en'or input signals are 
real, and are transmitted to the update generator 28. 

The update generator 28 includes junctions 54 and 60, multipliers 56, 58, 62 and 64. and summers 66 and 68. The 
set of m error input signals e from the in-phase element 50 of the en'or weighting element 26 is transmitted to junction 
54, where tiie signals e are split. From junction 54. one set of m error input signals e is provided to multiplier 56, and 
another set of m error input signals e is provided to multiplier 58. Likewise, the set of m error input signals e from the 
quadrature element 52 of the error weighting element 26 is transmitted to junction 60. where tiie signals e are split. From 
junction 60, one set of m en-or input signals e is provided to multiplier 62. and another set of m error input signals e is 
provided to multiplier 64. 

The m error input signals e provided to multiplier 62 are multiplied by the in-phase demodulation signal 70, which 
is preferably the same as the normalized in-phase reference signal 36. The m error input signals e provided to multiplier 
56 are multiplied by the quadrature demodulation signal 72, which is preferably tiie same as tiie normalized phase- 
shifted quadrature reference signal in line 43. This demodulation should occur during each sample period of adaptation. 
The output from nujltipliers 56 and 62 is summed in summer 66 to generate the negative of m updates Uj for the quadrature 
scaling vector Y| in the quadrature scaling element 44 that generates the quadrature reference signals yj. 



C = USV' 
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The m error input signals e provided to nrtuttiplier 58 are multiplied by the normalized in-phase demodulation signal 
76. The m error input signals e provided to multiplier 64 are multiplied by tiie normalized quadrature demodulation signal 
74. This demodulation should occur during each sample period of adaptation. The output from multipliers 58 and 64 is 
subtractively summed in summer 68 to generate m updates for the in-phase scaling vector Yp In the In-phase scaling 

5 element 40 that generates the m In-phase reference signals y^ 

The scaling vectors Yp and Y| are the adaptive parameters In the adaptive parameter bank 13. A reasonable bound 
on adaptation step size ^max is 0.25 divided by the number of sample periods corresponding to the average propagation 
delay tiirough the auxiliary path 24 between the actuators 16 and the error sensors 22. In cases where the error path 
24 is highly resonant, the step size \i should be smaller. 

10 The weighted en-or method of adaptation as shown in FIG. 2 approximates a real-time system and accounts for tiie 
propagation delay in the speaker-en-or patii 24 (I.e. auxiliary C path) because the updates are being accumulated over 
time. In the system 1 0 of FIG. 2. the in-phase and quadrature yj output signals are orthogonal and the update signals 
Ur and Uj from summers 66 and 68 provide adaptation essentially along these orthogonal directions. 

It is not necessary, however, that the con-ection signals y^ be generated from tiie combination of in-phase y^ and 

15 quadrature yj output signals as shown in FIG, 2. Referring to FIG. 3, the plurality of n (i.e. n = m ) correction signals yp 
can be generated by a parameter bank 1 3 in the form Acos(ci)k + <I>) where o) is the frequency estimate from the phase 
locked loop circuit 32. A is an adaptive amplitude, and O is an adaptive phase shift. In this case, the in-phase reference 
signals 70 and 72 are of the form cos(o)k + <P) , and the quadrature reference signals 72 and 74 are of the form 
sin(a)k + 4>) . The updates u^ from summer 68 update the amplitude A, and the ufbdates u^, from summer 66 update 

20 the phase shift <t> The amplitude A and the phase shift <t> are not in general orthogonal, such as the In-phase y^ and 
quadrature yj input signals shown in FIG. 2, and a system 10 as shown in FIG. 3 Is less likely to adapt along the shortest 
patii of the en-or surface as tiie system 1 0 shown in FIG. 2. 

The system 10 shown in FIG. 4 implements a constrained output control method. The system 10 in FIG. 4 uses an 
update generator 28 and an error weighting element 26 tiiat can be similar to the system 10 in FIG. 2, but the adaptive 

25 parameter bank 13 and the output weighting element 1 4 are preferably different than in FIG. 2. In Fig. 4, an input signal 
x(k) from an input sensor 30 is transmitted to the adaptive parameter bank 13 in the controller 12. The input signal x(k) 
is received by a phase locked loop circuit 32A, which Is now illustrated in FIG. 4 to Include a phase shifter such as phase 
shifter 42 shown in FIG. 2. The phase locked loop circuit 32A transmits an in-phase reference signal cos(cok) to junction 
80. Thereafter, the in-phase reference signal Is transmitted to scaling element 84, and to scaling element 92. Hie phase 

30 locked loop circuit 32A also transmits a quadrature reference signal sin(Q)k) to junction 82. Thereafter, tiie quadrature 
reference signal is transmitted to scaling element 90. and to scaling element 86. Scaling elements 84 and 90 use tiie 
same or a copy of the same adaptive in-phase scaling vector Yr. Likewise, scaling elements 86 and 92 use the same 
or a copy of the same adaptive quadrature scaling vector Y,, The adaptive parameters in the adaptive parameter bank 
13, thus includes an adaptive quadrature scaling vector Y|, an adaptive in-phase scaling vector Yr and copies of tiie 

35 same. 

The output from tiie scaling elements 84, 86. 90 and 92 is in general a set of 4 x m output signals y The output 
signals y are transmitted to the output weighting element 14 for processing In which the output signals are preferably 
scaled and combined to generate a set of n scaled and phase shifted correction signals y^- "I^e output signal from 
scaling element 84 Is subtractively summed with tiie output signal from ttie scaling element 86 in summer 88. The 
40 negative of the output signal from scaling element 90 is summed with the negative of the output signal from the scaling 
element 92 in summer 94. The output from tiie summer 88 is in general a set of m In-phase processing signals zp,e (i.e. 
real component). The output from the summer 94 is in general a set of m quadrature processing signals Zrq ("•©• imaginary 
component). 

The m in-phase processing signals zr^ from summer 88 are transmitted to an in-phase element 96 of the output 
45 weighting element 14, Re{Hi(ej'^)}. The m quadrature processing signals Zrq from summer 94 are transmitted to a 
quadrature element 98 of the output weighting element 1 4, lm{Hi (ei"^}. The in-phase element 96 of the output weighting 
element 14 contains the coefficients of the real parts of the complex elements of the output weighting matrix, Hi(ej'^). 
The quadrature element 98 of the output weighting element 14 contains the coefficients of the imaginary parts of tiie 
complex elements of the output weighting matrix, Hi (e*"). The in-phase 96 and the quadrature 98 elements of the output 
50 weighting element contain real values. Both the in-phase 96 and the quadrature 98 elements of the output weighting 
element 1 4 are preferably n x m matrices. The output from tiie in-phase 96 and the quadrature 98 elements of the output 
weighting element are summed in summer 1 00 to form n correction signals yp. The n conection signals yn are transmitted 
to an an-ay of n actuators 1 6 as discussed above. 

The structure of the adaptive parameter bank 13 and tiie output weighting element 14 represented as H^i^"^ can 
55 be understood by expanding the term H^y into complex components (i.e., 
H 1 y = (Re {H 1 } y , - Im {H 1 } y ,) + i(lm {H 1 ) y , + Re {H 1 } y ,) ). It is preferred that (e*") depend on the C model to 
facilitate convergence of the system to a stable solution. The output weighting element 14 therefore scales, phase-shifts 
and combines output signals y to cause the system to converge and to selectively constrain tiie correction signals yn- 
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The remaining part of the system 1 0 shown in FIG. 4 can be similar to the system 1 0 desaibed in FIG. 2 except for 
certain distinctions. For instance, the in-phase update signal Up from summer 68 is used to adapt the scaling elements 
84 and 90 (i.e., adapt vector YfO- The quadrature update signal Uj from summer 66 is used to adapt the scaling elements 
86 and 92 (i.e., adapt vector Y|). 

5 In FIG. 4. the output from the controller 12 to the n actuators (i.e. the n correction signals yj can be constrained 
using the techniques of singular value decomposition as described above. Preferably, this can be done by setting the 
output weighting element 14 (i.e. output weighting matrix Hi(e'^)) to the n x n matrix V. and the error weighting element 
26 (i.e. error weighting matrix H2(e'^)) to -N^U^ where is the Hermitian transpose of the p x p U matrix, and is 
the Hermitian transpose of the p x n normalizing matrix N which is formed by taking the transpose of S and inverting 

10 some of the values along the diagonal S (e.g., the values that are not zero or close to zero). In such a system, the 
secondary input does not contain components corresponding to the effective null space of C. That is. the secondary 
input does not contain components corresponding to the columns in V corresponding to singular values in C. Processing 
selected output signals yj and yr corresponding to values in S that are not too close to zero through V to generate the 
correction signals thus eliminates undesirable components in y^ that may cause problems associated with over param- 

15 eterization. Also, processing the error signals ep through -N^U^ restricts adaptation corresponding the null elements of 
N (i.e. the effective null space of C). 

If none of the diagonal elements in N are set to zero, then in the case where there are more en'or sensors 22 than 
actuators 16, i.e. p > n, the system converges to the same minimum cost function where the solution is given by: 

20 ynM-<c''cr'c''d (3) 

where is the Hermitian transpose of the p x n C matrix, and d is the disturbance designated by reference numeral 
18. In addition, the system 10 converges in less time than a system implementing the gradient descent method. In the 
case where tiiere are less error sensors 22 than actuators 1 6, i.e. p < n, tiie controller output yn represents tiie minimum 

25 required energy in the n secondary inputs for complete cancellation of the disturbance 1 8. 

The output 1 4 and the en'or 26 weighting elements can be chosen to rotate the coordinates of the system such that 
each of the n processing signal pairs (i.e., n pairs of Zp^ and Zpe) affects adaptation along one of the principle axis of 
the hyper-elliptical quadratic surfaces of the cost function for the correction signals y^. These principal axes are ortiiog- 
onal. Adaption along these axes depends on the steepness of the axis; however, normalizing with matrix N normalizes 

30 adaptation with respect each axis such that adaptation proceeds along each axis at the same rate. Therefore, adaptation 
occurs along a straight line in the controller parameter space towards an optimum solution. Further, the system is decou- 
pled in that each of the n pairs of Zp^ and Zp^ do not affect the adaptation of the otiier pairs of Zp^ and Zp^. 

Moreover, using the output weighting element 14 to provide constrained output allows for improved control over 
parameter leakage and for improved control over system nrx)nitoring because it allows selective leakage or monitoring 

35 along selected axes. In addition, the controller 1 2 can be constrained from adapting along axes corresponding to singular 
values by testing for small singular values and setting the corresponding elements of mati'ix N to zero, or equivalentiy, 
by not adapting the con^esponding processing signals Zp^ or Zpg. 

When tracking, it may be useful to switch between an old output weighting element 14 (i.e. old output weighting 
matrix Hi(e'^oy) to a new output weighting element 14 (i.e. new output weighting matrix HT(e''^new)- the preferred 

40 case where H ^(e^'^) = V , the transition from an old set of correction signals y^ oid to a new set of correction signals 
yn.new be Smoothed by using V^pew Vqw yn.ow as the first yn,nQvt this manner, the set of correction signals from 
the previous sample period yn,oid "S matched to the set of correct signals for the present sample period yn.new 

Anotiier method to constrain the output from the controller 12 (i.e. constraining the n correction signals y^), using 
the techniques of singular value decomposition, is to omit the error weighting element 26 (i.e. set to identity), and set 

45 the output weighting element 14 to -VN^U" as shown in FIG. 5. In tiie case of n = p where none of the values of n are 
nulled (i.e. all of the singular values are sufficientiy large). H^ = - VN^u" is equal to •C""'. An advantage of this method 
is that the output from tiie controller 1 2 can be consti-ained, and only a single processing matrix is required. Alternatively, 
the output weighting element 1 4 can apply an n x p matrix H ^ = -C " where is the Hermitian transpose of the p x n 
C matrix. 

50 It is recognized that various equivalents, alternatives and modifications are possible within the scope of the appended 
claims. 

Claims 

55 1 . An adaptive tonal control system having a system input and a system output, the adaptive tonal control system 
comprising: 

a plurality of actuators each receiving a correction signal and outputling a secondary input, the secondary 
inputs combining with the system input to yield the system output; 

a plurality of error sensors sensing the system output, each error sensor generating an error signal that can 
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be used to update adaptive parameters: and 

an adaptive controller that outputs the correction signals, the controller having 

an adaptive parameter bank that outputs a plurality of output signals in accordance with adaptive parameters, 
a C model of a path between the output of the adaptive controller and the error sensors, the C model having 

an effective range and an effective null space, and 

an output weighting element that inputs the output signals from the parameter bank and weights the output 

signals to generate correction signals which are constrained to be wrthin the effective range of the C model. 

2. A system as recited in claim 1 wherein the effective range of the C model does not include the effective null space 
of the C model, and the effective null space of the C model is defined by iCyj being small, where C is a matrix 
representing the C model, and yn is a vector representing correction signals having non-trivial values. 

3. A system as recited in claim 1 wherein: the C model can be represented as a p x n matrix C that can be decomposed 
as: C = USV " , where S is a p x n matrix in which the off diagonal elements are zero. U is a p x p unitary matrix, 
and is the Hermitian transpose of an n x n unitary matrix V; and 

the output weighting element can include a matrix representing -VS^U", where is the Hermitian transpose 
of S and is the Hermitian transpose of U. 

4. A system as recited in claim 3 wherein one or more of the diagonal elements in are normalized. 

5. A system as recited in daim 1 further comprising an error weighting element that inputs the error signals from the 
error sensors and weights the en-or signals to generate error input signals. 

6. A system as recited in claim 5 wherein: 

the C model can be represented as a p x n matrix C that can be decomposed as: C = USV " , where S is a 
p x n matrix in which the off diagonal elements are zero, U is a p x p unitary matrix, and is the Hermitian transpose 
of an n X n unitary matrix V; 

the output weighting element includes a matrix representing V; and 

the error weighting element includes a matrix representing -S^U^, where S*^ is the Hermitian transpose of 
S and is the Hermitian transpose of U. 

7. A system as recited in claim 6 wherein one or more of the diagonal elements in are normalized. 

8. A system as recited in claim 1 wherein the adaptive parameter bank outputs a plurality of in-phase output signals, 
and a plurality of quadrature output signals. 

9. A system as recited in claim 1 further comprising an input sensor that senses the system input and generates a 
reference signal in response thereto. 

1 0. A system as recited in claim 9 further comprising a phase locked loop circuit that receives the reference signal from 
the input sensor and outputs an in-phase reference signal that is used in the adaptive parameter bank. 

1 1. A system as recited in claim 10 further comprising a phase shifter that receives a copy of the in-phase reference 
signal and outputs a quadrature reference signal that is 90** out-of-phase from the in-phase reference signal and is 
also used in the adaptive parameter bank. 

12. A system as recited in claim 1 1 wherein the output weighting element has an in-phase output weighting element, 
and a quadrature output weighting element, and an output summer that sums a signal from the in-phase output 
weighting element and from the quadrature output weighting element and outputs the correction signals. 

1 3. An adaptive tonal control system having a system input and a system output, the adaptive control system comprising: 

a plurality of actuators each receiving a correction signal and outputting a secondary input, the secondary 
inputs combining with the system input to yield the system output; 

a plurality of error sensors sensing the system output, each sensor generating an error signal; 
an adaptive controller that outputs the correction signals, the controller having 

an adaptive parameter bank that outputs a plurality of output signals in accordance with the adaptive param- 
eters, the output signals being used to generate the correction signals, 

a C model of a path between the output of the adaptive controller and the error sensors, the C model having 
an effective range and an effective null space. 
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an error weighting element that Inputs the enor signals from the error sensors and weights the error signals 
to generate en-or input signals that are used to update the adaptive parameters in the adaptive parameter bank 
such that adaptation is constrained to be within the effective range of the C model. 

5 1 4. A system as recited in daim 1 3 wherein the effective range of the C model does not include the effective null space 
of the C model and the effective null space of the C model is defined by iCynQ being small, where C is a matrix 
representing the C model, and yn is a vector representing con'ection signal having non-trivial values. 

1 5. A system as recited in daim 1 3 wherein the C model can be represented as a p x n matrix C that can be decomposed 
10 as: C - USV " . where S is a p x n matrix in which the off diagonal elements are zero, U is a p x p unitary matrix, 
and is the Hermitian transpose of an n x n unitary matrix V; and 

the error weighting element can include a matrix representing -VS^U^. where S*^ is the Hermitian transpose 
of S and is the Hermitian transpose of U. 

15 1 6. A system as recited in claim 1 5 wherein one or more of the diagonal elements in are normalized. 

17. A system as recKed in claim 13 wherein each correction signal is the summation of an in-phase output signal and 
a quadrature output signal. 

20 18. A system as recited in claim 13 further comprising an input sensor that senses the system input and generates a 
reference signal. 

1 9. A system as recited in daim 1 8 further comprising a phase locked loop drcuit that receives the signal from the input 
sensor and outputs an in-phase reference signal. 

25 

20. A system as recited in claim 19 further comprising a phase shifter that can receive a copy of the in-phase reference 
signal and output a quadrature reference signal. 

21 . A system as recited in daim 1 3 wherein the error weighting element has an in-phase element that receives the error 
30 signals and outputs an in-phase error input signals, and a quadrature element that receives the error signals and 

outputs quadrature error input signals. 

22. A system as recited in claim 21 further comprising a parameter update generator that applies an in-phase demod- 
ulation signal and a quadrature demodulation signal to the in-phase and quadrature error input signals every sample 

35 period, and outputs in-phase update signals and quadrature update signals which are used to adapt the adaptive 
parameters in the adaptive parameter bank. 

23. An adaptive tonal control system having a system input and a system output, the adaptive control system comprising: 

a plurality of actuators, each receiving a con-ection signal and outputting a secondary input, the secondary 
40 inputs combining with the system input to yield the system output: 

a plurality of error sensors sensing the system output, each sensor generating an en-a signal; and 
an adaptive controller that outputs the correction signals and has 

an adaptive parameter bank that inputs a reference signal during sequential sample periods, and outputs a 
plurality of output signals in accordance with adaptive parameters for each of the sequential sample periods, the 
45 output signals being used to generate the con-ection signals, 

an error weighting element that receives the error signals and outputs a plurality of error input signals, 
a parameter update generator that applies an in-phase demodulation signal and a quadrature demodulation 
signal to the error input signals for each sample period when generating update signals, and outputs an in-phase 
update signal and a quadrature update signal that are used to update the adaptive parameters in the adaptive 
50 parameter bank. 

24. A system as recited in claim 23 wherein: 

the error-weighting element has an in-phase weighting element that receives the enror signals and outputs 
in-phase error input signals and a quadrature element that receives the en-or signals and outputs quadrature error 
55 input signals; and 

the parameter update generator has a first multiplier that receives quadrature error input signals and the in- 
phase demodulation signal and outputs a first signal to a quadrature summer, a second multiplier that receives the 
in-phase error input signals and the quadrature demodulation signal and outputs a second signal to the quadrature 
summer, the quadrature summer outputting the quadrature update signal, a third multiplier that receives the quad- 
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rature error input signals and the quadrature demodulation signal and outputs a third signal to an ini3hase summer, 
and a fourth multiplier that receives the in-phase error input signals and the in-phase demodulation signal and 
outputs a fourth signal to the in-phase summer, the in-phase summer outputting the in-phase update signal. 

25. In an adaptive control system having a system input and a system output, the method of controlling a tone comprising 
the steps of: 

modeling a path between an output of an adaptive filter controller and one or more error sensors in the system 
to generate a C path model; 

generating a plurality of output signals: 

weighting the plurality of output signals to yield a plurality of correction signals that are within the effective 
range of the C path model; and 

using the correction signals to provide a secondary input that combines with the system input to yield the 
system output. 

26. A method as recited in claim 25 wherein the effective range of the C path model is determined by nulling components 
of the correction signals in which gCynO is small, where C is a matrix representing the C path model and yn is a vector 
representing non-trivial correction signals. 

27. A method as recited in daim 25 further comprising the steps of: 

sensing the system output with the error sensors to generate a plurality of error signals; 

using the error signals to update the adaptive parameters that are used to generate the output signals. 

28. A method as recited in daim 27 further comprising the step of: 

weighting the plurality of error signals to yield a plurality of error input signals; and 
using the error input signals to update the adaptive parameters. 

29. A method as recited in daim 27 further comprising the steps of: 

tracking the frequency of the tone being controlled; 

matching the plurality of correction signals of a previous sample periods to the plurality of correction signals 
for a present sample period to determine the weighting of the correction signals of the present sample period. 

30. In an adaptive control system having a system input and a system output, a method of controlling a tone conprising 
the steps of: 

modeling a path between an output of an adaptive controller and one or more error sensors in the system to 
generate a C path model; 

generating a plurality of output signals in accordance with adaptive parameters; 

using the output signals to provide a secondary input to the system which combines with the system input 
to yield the system output; 

sensing the system output with error sensors to generate a plurality of error signals; 

weighting the plurality of error signals to yield a plurality of error input signals such that adaptation is con- 
strained with respect to the effective range of the C model; and 

using the error input signals to update the adaptive parameters. 

31 . A method as recrted in daim 30 wherein the effective range of the C path model is determined by nulling components 
of the correction signals in which iCyj is small, where C is a matrix representing the C path model and y^ is a vector 
representing non-trivial correction signals. 

32. A method as recited in daim 30 wherein: 

the plurality of error signals are weighted to yield the plurality of in-phase error input signals and a plurality 
of quadrature error input signals; and 

the adaptive parameters are updated using an in-phase update signal and a quadrature update signal that 
are generated by modulating the in-phase en-or input signal and the quadrature error input signal during the sample 
periods when adaptation is occurring. 

33. A method as recited in daim 30 further comprising the step of: 

weighting the plurality of output signals to yield a plurality of correction signals that are within the effective 
range of the C path model, and using the correction signals to provide the secondary input that combines with the 
system input to yield the system output. 
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34. A method as recited in claim 30 further comprising the step of leaking components of the output signals in the 
effective null space of the C model. 

35. In an adaptive control system having a system input and a system output, a method of controlling a tone comprising 
5 to steps of: 

inputting a reference signal for sequential sample periods to an adaptive controller: 

outputting a plurality of correction signals from the adaptive controller which are used to generate secondary 

inputs that combine with the system input to yield the system output; 

sensing the system output with error sensors to generate a plurality of error signals; 
10 modeling a path between the output of the adaptive controller and the error sensors to generate a C path 

model; 

weighting the error signals to generate a plurality of error input signals; 

applying an in-phase demodulation signal and a quadrature demodulation signal to the en-or input signals 
every sannple period to generate an In-phase update signal and a quadrature update signal for sample periods when 
15 the system is generating update signals; 

using the in-phase and the quadrature update signals to update the adaptive parameters. 

36. The method as recited in claim 35 further comprising sensing the system input with input sensors to generate one 
or more reference signals every sample period. 

20 

37. A method as recited in claim 36 further comprising the step of processing each reference signal to generate an in- 
phase reference signal and a quadrature reference signal. 

38. An adaptive tonal control system having a system input and a system output, the adaptive tonal control system 
25 comprising: 

a plurality of actuators each receiving a correction signal and outputting a secondary input, the secondary 

inputs combining with the system input to yield the system output; 

a plurality of en-or sensors sensing the system output, each error sensor generating an error signal that can 

be used to update adaptive parameters; and 
30 an adaptive controller that outputs the correction signals, the controller having 

an adaptive parameter bank that outputs a plurality of output signals in accordance with adaptive parameters, 

a C model of a path between the output of the adaptive controller and the error sensors, and 

an output weighting element that inputs the output signals from the parameter bank and weights the output 

signals to generate correction signals, wherein the output weighting element depends on the C model. 

35 

39. A system as recited in claim 38 further comprising an error weighting element that inputs the error signals from the 
error sensors and weights the error signals to generate error input signals. 

40. A system as recited in claim 38 wherein the adaptive parameter bank outputs a plurality of in-phase output signals, 
40 and a plurality of quadrature output signals. 

41 . A system as recited in claim 38 further comprising an input sensor that senses the system input and generates a 
reference signal in response thereto. 

45 42. A system as recited in daim 41 further comprising a phase locked loop circuit that receives the reference signal 
from the input sensor and outputs an in-phase reference signal that is used in the adaptive parameter bank. 

43. A system as recited in claim 42 further comprising a phase shifter that receives a copy of the in-phase reference 
signal and outputs a quadrature reference signal that is 90" out-of -phase f rom the in-phase reference signal and is 

50 also used in the adaptive parameter bank. 

44. A system as recited in claim 43 wherein the output weighting element has an In-phase output weighting element 
that receives an in-phase output signal from the adaptive parameter bank, and a quadrature output weighting ele- 
ment that receives a quadrature output signal from the adaptive parameter bank, and an output summer that sums 

55 a signal from the in-phase output weighting element and from the quadrature output weighting element and outputs 
the correction signals. 

45. In an adaptive control system having a system input and a system output, the method of controlling a tone comprising 
the steps of: 
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modeling a path between an output of an adaptive filter cx)ntrolter and one or more error sensors in the system 
to generate a C path model; 

generating a plurality of output signals; 

processing the plurality of output signals to yield a plurality of correction signals In such a manner that the 
processing depends on the C model so that the system converges to a stable solution; and 

using the correction signals to provide a secondary Input that combines wrth the system input to yield the 
system output. 

46. A method as recited in daim 45 wherein the processing is accomplishing by filtering. 

47. A method as recited in daim 45 wherein the processing of the output signals indudes scaling, phase-shifting and 
combining. 

48. A method as recited in daim 45 wherein the output signals are generated by; 

processing an in-phase reference signal with an in-phase scaling vector Yp to generate a first plurality of 
signals; 

processing the In-phase reference signal with a quadrature scaling vector Y| to generate a second plurality 
of output signals; 

processing a quadrature reference signal with the In-phase scaling vector to generate a third plurality of 
output signals; and 

processing the quadrature reference signal with the quadrature scaling vector Y) to generate a fourth plurality 
of output signals. 

49. A method as recited in claim 46 wherein the four output signals are processed to generate the plurality of correction 
signals by; 

summing the first and the fourth plurality of output signals to form a plurality of in-phase processing signals 

summing the second and the third plurality of output signals to form a plurality of quadrature processing 
signals ztm; 

processing the in-phase processing signals zr^ through an in-phase output weighting element that depends 
on the C model; 

processing the quadrature processing signals Zif^, through a quadrature output weighting element that 
depends on the C model; and 

summing the plurality of processed in-phase processing signals and the plurality of processed quadrature 
processing signals to form the plurality of correction signals. 

50. A method of controlling tones in an acoustic system, the method comprising providing an adaptive controller respon- 
sive to an error signal representing system output for providing a correction signal to be combined with a system 
input, the adaptive controller including an adaptive model of a path between the output of the controller and the 
system output, characterised In that the method includes the step of constraining the controller output with respect 
to the null space (or effective null space) of the model. 
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